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Abstract 

We consider QGD radiative corrections to the prodnction of fonr charged leptons in 
the ZZ signal region at the LHG. We report on the complete calcnlation of the next-to- 
next-to-leading order (NNLO) corrections to this process in QGD pertnrbation theory. 
Nnmerical resnlts are presented for ^fs = 8 TeV, nsing typical selection cnts applied 
by the ATLAS and GMS collaborations. The NNLO corrections increase the NLO 
hdncial cross section by abont 15%, and they have a relatively small impact on the 
shape of the considered kinematical distribntions. In the case of the A$ distribntion of 
the two Z candidates, the NNLO corrections improve the agreement of the theoretical 
prediction with the GMS data. 
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The production of Z-boson pairs at the Large Hadron Collider (LHC) provides an important 
test of the electroweak (EW) sector of the Standard Model (SM) at the TeV scale. Small deviations 
in the observed rates or in the kinematical distributions could be a signal of new physics, possibly in 
the form of anomalous couplings. At the same time, ZZ production is an irreducible background 
for Higgs boson production and new-physics searches. Particularly important are the off-shell 
effects below the ZZ threshold, relevant for the Higgs signal region in the four-lepton channel. 
Various measurements of ZZ hadropro duct ion have been carried out at the Tevatron and the LHC 
(for some recent results see Refs. [IHZ])- 

From the theory side the hrst NLO predictions for on-shell ZZ production were obtained long 
ago [Slin] • The leptonic decays of the Z bosons were included, initially neglecting spin correlations 
in the virtual contributions HD]. The computation of the relevant one-loop helicity amplitudes m 
enabled the hrst complete NLO calculations DM, including spin correlations and off-shell effects. 
The loop-induced gluon-fusion production channel, which formally contributes only at the next- 
to-next-to-leading order (NNLO), was computed in Refs. [TTlfTS] . The corresponding leptonic 
decays were taken into account in Refs. [TGHIH]. NLO predictions for ZZ production including 
the gluon-induced contribution, the leptonic decays with spin correlations and off-shell effects were 
presented in Ref. [I9]. The NLO QCD corrections to on-shell ZZ-|-jet production were discussed 
in Ref. [2011^ . and the EW corrections to ZZ production in Ref. [22]. A decisive step forward was 
carried out in Ref. [23] where the inclusive NNLO cross section for on-shell ZZ production was 
presented. This calculation was based on the evaluation of the two-loop amplitude for on-shell 
ZZ production. Later, the two-loop helicity amplitudes for all the vector-boson pair production 
processes were presented [211I2S]- This computation paves the way to the consistent inclusion of 
the leptonic decays and off-shell effects in the NNLO computation. 

In this Letter we carry out this step by considering ZZ production at NNLO including the 
leptonic decays of the vector bosons together with spin correlations and off-shell effects. Contri¬ 
butions from Zy* and 7 * 7 * production as well as from pp —)■ Z/ 7 * —)■ 4 leptons topologies are also 
consistently included with all interference terms. Our calculation allows us to apply arbitrary cuts 
on the hnal-state leptons and the associated QCD radiation. Here we present selected numerical 
results for pp —)■ 4 leptons at the LHC in NNLO QCD, using the typical cuts that are applied in 
the experimental ZZ analyses. 

Our calculation is performed with the numerical program MatrixQ, which combines the qt 
subtraction [26] and resummation [2^ formalisms with the Munich Monte Carlo framework [28] . 
Munich provides a fully automated implementation of the Catani-Seymour dipole subtraction 
method [23130], an efficient phase-space integration, as well as an interface to the one-loop gener¬ 
ator OpenLoops [31] to obtain all required (spin- and colour-correlated) tree-level and one-loop 
amplitudes. For the numerically stable evaluation of tensor integrals, OpenLoops relies on the 
Collier library [32], which is based on the Denner-Dittmaier reduction techniques [331134] and the 
scalar integrals of [33]. To deal with problematic phase-space points, a rescue system is provided, 
which employs the quadruple-precision implementation of the OPP method in CutTools [3S] 
and scalar integrals from OneLOop [37]. Our implementation of qt subtraction and resummation 
for the production of colourless hnal states is fully general, and it is based on the universality of 
the hard-collinear coefficients [38] appearing in transverse-momentum resummation. These coef- 

'I' Matrix is the abbreviation of “Munich Automates qT subtraction and Resummation to Integrate Cross 
Sections”, by M. Grazzini, S. Kallweit, D. Rathlev, M. Wiesemann. In preparation. 
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ficients were explicitly computed for quark-initiated processes in Refs. [3SH1I]- For the two-loop 
helicity amplitudes we use the results of Ref. |25], and of Ref. |12] for Drell-Yan like topologies. 

A preliminary version of MATRIX has been employed in the NNLO computations of Refs. [231 
HSHIS] , and in the resummed calculation of Ref. |16] . 

We consider pp collisions with y/s = 8 TeV. As for the EW couplings, we use the so-called 
Gfj, scheme, where the input parameters are Gp, mw, ^z- More precisely, consistent with the 
OpenLoops implementation, we use the complex W and Z boson masses to dehne the EW 
mixing angle as cos^^ = (m^ — iVw'^w)lijnz ~ i^z'^nz)- In particular, we use the values 
Gf = 1.16639 X 10-5 GeV-2, mw = 80.399 GeV, Tw = 2.1054 GeV, mz = 91.1876 GeV, 
Tz = 2.4952 GeV. For the top quark we use mt = 173.2 GeV, Ft = 1.4426 GeV, and for the Higgs 
boson rnn = 125 GeV, Tp = 4.07 MeV. Both the top quark and the Higgs boson only appear in 
diagrams with closed top-quark loops, thus entering the gluon-fusion channel and the real-virtual 
contributionjil We use the NNPDF3.0 m sets of parton distributions with as{mz) = 0.118, and 
the as running is evaluated at each corresponding order (i.e., we use (n -|- l)-loop as at N”LO, 
with n = 0, 1,2). We consider Nf = 5 massless quark flavours. The central renormalization {pr) 
and factorization (/i^) scales are set to Pr = Pf = f^z- 

We hrst consider the ATLAS analysis of Ref. |5] in the three decay channels e’''e-e’''e-, 
and e~^e~p^p~. The invariant masses of the two reconstructed lepton pairs are 
required to fullhl the condition 66 GeV < mu <116 GeV. In the case of two lepton pairs with the 
same flavours there is a pairing ambiguity, which is resolved by choosing the pairing that makes 
the sum of the absolute distances from mz smaller. The leptons are required to have Pt > 7 GeV 
and rapidity \p\ < 2.7. For any lepton pair we require AR{l,l') > 0.2, independently of the 
flavours and charges of I and I'. 

The corresponding cross sections are reported in Tab. [H where the ATLAS results are also 
shown. The uncertainties on our theoretical predictions are obtained by varying the renor¬ 
malization and factorization scales in the range 0.5mz < Pr, Pf < 2 mz with the constraint 
0.5 < /ip/pi? < 2. Independently of the leptonic decay channels, the NNLO corrections increase 
the NLO result by about 15%, similarly to what was found for the inclusive cross section for 
on-shell ZZ production [23]. This is as expected because the selection cuts are mild and do not 
signihcantly alter the impact of radiative corrections. The scale uncertainties are about ±3% at 
NLO and remain of the same order at NNLO. As noted for the inclusive cross section [23], the 
NLO scale uncertainty does not cover the NNLO effect. This is not surprizing since the loop- 
induced gluon-fusion contribution, which provides ~ 60% of the O (ag) correction, opens up only 
at NNLO. The NNLO corrections improve the agreement of the theoretical prediction with the 
data for the channel, whereas they deteriorate the agreement in the case of the 4e and 

4/i channels. We note, however, that the predicted hducial cross sections are still consistent with 
the ATLAS measurements at the la level within the statistics-dominated uncertainties. 

Secondly, we consider the GMS analysis of Ref. [7]. The hducial region is dehned as follows: 
all muons are required to fulhll > 5 GeV, |p^| < 2.4, while all electrons are required to fulhll 
pfn > 7 GeV, |p®| < 2.5. In addition, the leading- and subleading-lepton transverse momenta must 
satisfy ph > 20 GeV and ph > 10 GeV, respectively. In the case of two lepton pairs with the 

^ The Higgs boson contributes less than 1% to the loop-induced gg —>■ ZZ cross section, whereas its effect on 
the real-virtual contribution is numerically negligible. 
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Ghannel 

ctlo (fb) 

O'NLO (fb) 

O'NNLO (fb) 

^exp (fb) 

e'’'e“e'’'e“ 

3.547(1)1^:^^ 

5.047(1)1^5 


4.65:?(stat)55syst.)5;}(lumi.) 

jjL'^ IJi~ fi,'^ fJi~ 

5.o5:5(stat)5:2(syst.)5;2(iumi-) 


6.950(1)75 

9.864(2)55 

11.31(2)55 

ll.l55stat)5;5(syst.)5;3(lumi.) 


Table 1: Fiducial cross sections and scale uncertainties for ATLAS cuts at LO, NLO, and NNLO 
in the three considered leptonic decay channels. The ATLAS data are also shown. 


same flavours, the pairing ambiguity is resolved by choosing the pair with the smallest distance 
from rriz- This pair is called Zi, the remaining pair is called Z2. The invariant masses of the two 
reconstructed lepton pairs are required to fulhll 60 GeV < mu < 120 GeV. We note that in the 
case of identical flavours this dehnition of the hducial region does not prevent the invariant masses 
of the other two possible lepton pairs from becoming arbitrarily small, giving rise to a collinear 
7* —)■ singularity. To avoid that, we follow GMS and add an additional cut mu > 4 GeV 

on all lepton pairs of the same flavours and opposite charges]! The corresponding hducial cross 
sections and scale uncertainties, computed as above, are reported in Tab. [2l Like for the ATLAS 
analysis, the NNLO corrections increase the NLO hducial cross section by about 15%. The scale 
uncertainties are similar to those reported in Tab. [U 


Ghannel 

ctlo (fb) 

cnlo (fb) 

c’'nnlo (fb) 

e'’'e“e'’'e“ 

3.149(1)75: 

4.493(1)75: 

5.16(1)5S 

/i+/r-/i+/r- 

2.973(1)7;« 

4.255(1)55 

4-90(1)5:^: 


6-179(l)lS 

8.822(1)75: 

10.15(2)75 


Table 2: Fiducial cross sections and scale uncertainties for GMS cuts at LO, NLO, and NNLO in 
the three considered leptonic decay channels. 

GMS does not report the hducial cross sections corresponding to the above cuts, but only 
normalized distributions, to which we compare our results. We start with the invariant-mass 
distribution of the four leptons, which is depicted in Fig. [H The lower panels show the theory/data 
comparison, and the NNLO result normalized to the central NLO prediction. We see that the 
NNLO corrections have a limited impact in the comparison with the data, which still have large 
uncertainties. The NNLO ehects on the normalized distribution are relatively small: they are 
completely negligible at low invariant masses, and they increase to —5% in the high mass region. 
This means that the NNLO corrections make the invariant mass distribution slightly softer. We 
have checked that this effect is due to the gluon-fusion contribution, whose relative effect decreases 
at high masses, due to the larger values of Bjorken x that are probed. The NLO (NNLO) scale 
uncertainties range from about ±2% (±1%) at low mzz to ±4% (±2%) at high mzz- 

^We thank Alexander Savin for providing us with this information. 
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Figure 1: The four-lepton invariant-mass distribution at NLO and NNLO compared to the CMS 
data. In the lower panels the ratio of our theoretical results over the data, and the NNLO result 
normalized to the central NLO prediction are presented. The bands correspond to scale variations 
as described in the text. 


4 
























































































Figure 2: The leading-lepton px (left) and the A0 (right) distributions at NLO and NNLO 
compared to the CMS data. In the lower panels the ratio of our theoretical results over the 
data, and the NNLO result normalized to the central NLO prediction are presented. The bands 
correspond to scale variations as described in the text. 


In Fig. |2]we show the analogous results for the leading-lepton px distribution (left) and the 
azimuthal separation (Ad)) of the two Z candidates (right). As in Fig. [1], we see that the NNLO 
effects on the px distribution do not change the comparison with the data in a signihcant way. 
The NNLO corrections are also relatively small in most of the range of px considered, except 
for the low px region, where they increase significantly. This effect is due to the gluon-fusion 
contribution, whose relative impact increases as px decreases. The situation is different for the 
Ad) distribution. Here the NNLO corrections improve the agreement with the data, except for 
the first bin, where the CMS measurement is an order of magnitude below the theoretical NNLO 
prediction. The larger impact of NNLO corrections in the Ad) distribution can be understood 
easily by the observation that at LO the reconstructed Z bosons are always back-to-back, i.e., 
Ad)(Zi,Z2) = TT. As a consequence, the NNLO calculation is effectively NLO in the region 
0 < Ad) < TT. The NNLO corrections amount to about -|-25% when Ad) <1.5, and decrease as 
Ad) increases. We note that this effect is entirely due to the NNLO corrections to the qq channel 
addressed in this paper, since the loop-induced gluon-fusion contribution, which also enters at 
NNLO, affects the Ad) distribution only at Ad) = vr. The NLO scale uncertainties are about 
±11%, while at NNLO the uncertainties are about ±5% at low Ad), and decrease to about ±2% 
at high Ad). 

We have presented the hrst complete NNLO QCD calculation for the production of four charged 
leptons in the ZZ signal region at the LHC. We have studied the impact of NNLO corrections on 
the hducial cross sections and distributions measured by ATLAS and CMS at the LHC. As for 
the hducial cross sections, we found about ±15% NNLO corrections w.r.t. the NLO prediction, 
consistent with what was found for the inclusive cross section for on-shell ZZ production j23j . 
The impact on the normalized distributions we considered is small compared to the experimental 
uncertainties, but leads to an improved agreement with the data in the case of the Ad) distribution 
of the two Z candidates. Our calculation was performed with the numerical program Matrix, 
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which is able to carry out fully exclusive NNLO computations for a wide class of processes at 
hadron colliders. We look forward to further applications of our framework to other important 
LHC processes. 
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